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Binary pulsars provide an excellent system for testing general relativity because of their intrinsic rota- 
tional stability and the precision with which radio observations can be used to determine their orbital 
dynamics. Measurements of the rate of orbital decay of two pulsars have been shown to be consistent 
with the emission of gravitational waves as predicted by general relativity 1 2 , providing the most con- 
vincing evidence for the self-consistency of the theory to date. However, independent verification of 
the orbital geometry in these systems was not possible. Such verification may be obtained by deter- 
mining the orientation of a binary pulsar system using only classical geometric constraints, permitting 
an independent prediction of general relativistic effects. Here we report high-precision timing of the 
nearby binary millisecond pulsar PSR J0437— 4715, which establish the three-dimensional structure 
of its orbit. We see the expected retardation of the pulse signal arising from the curvature of space- 
time in the vicinity of the companion object (the 'Shapiro delay'), and we determine the mass of the 
pulsar and its white dwarf companion. Such mass determinations contribute to our understanding of 
the origin and evolution of neutron stars 3 . 



Discovered in the Parkes 70-cm survey 4 , 
PSR J0437— 4715 remains the closest and brightest 
millisecond pulsar known. It is bound to a low- 
mass helium white dwarf companion 5 ' 6 in a nearly 
circular orbit. Owing to its proximity, relative mo- 
tion between the binary system and the Earth sig- 
nificantly alters the line-of-sight direction to the pul- 
sar and, consequently, the orientation of the basis 
vectors used in the timing model (see Fig. 1). Al- 
though the physical orientation of the orbit in space 
remains constant, its parameters are measured with 
respect to this time-dependent basis and therefore 
also vary with time. Variations of the inclination an- 
gle, i, change the projection of the semi-major axis 
along the line-of-sight, x = a p sini/c, where a p is the 
semi-major axis of the pulsar orbit. 

The heliocentric motion of the Earth induces a 
periodic variation of x known as the annual-orbital 
parallax 7 , 

x° hs (t)=x int [l + ^r (B (t)-n'}. (1) 
The superscripts 'obs' and 'int' refer to the observed 



and intrinsic values, respectively, r®(i) is the posi- 
tion vector of the Earth with respect to the barycen- 
tre of the Solar System as a function of time, d is the 
distance to the pulsar, and tt' = sinfJIn — cosfiJo 
(see Fig. 1). Similarly, the proper motion of the bi- 
nary system induces secular evolution of the pro- 
jected semi-major axis 8 9 , such that: 

x ohs = x Ult -xcoti n -Sl'; (2) 

where /x = /i Q Io + M<5 Jo is the proper motion vector 
with components in right ascension, \i a , and declina- 
tion, fig. An apparent transverse quadratic Doppler 
effect (known as the Shklovskii effect) also arises 
from the system's proper motion and contributes to 
the observed orbital period derivative 10 : 

P£ hs = Pi nt + pP h ; (3) 

where /3 = p?djc, and [i = 

Observations of PSR J0437— 4715 were con- 
ducted from 11 July 1997 to 13 December 2000, us- 
ing the Parkes 64 m radio telescope. Over 50 ter- 
abytes of baseband data have been recorded with 
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the S2 Recorder 11 and the Caltech Parkes Swinburne 
Recorder (CPSR) 12 , followed by offline reduction 
at Swinburne's supercomputing facilities. Average 
pulse profiles from hour-long integrations were fit- 
ted to a high signal-to-noise template 13 , producing 
a total of 617 pulse arrival time measurements with 
estimated errors on the order of 100 ns. 

Previously considered negligible, the annual- 
orbital parallax has been largely ignored in exper- 
imental time-of-arrival analyses to date. However, 
our initial estimates of its peak-to-peak amplitude 
for PSR J0437-4715 (~ 400 ns) demonstrated that 
it would be clearly detectable above the timing 
noise. As can be seen in equation a; obs varies 
with a period of one year and phase determined by 
fl'. Its inclusion in our timing model therefore pro- 
vides a geometric constraint on fi. We also note 
that the value of i obs = (7.88 ± 0.01) x 10~ 14 ob- 
served in our preliminary studies is many orders 
of magnitude larger than the intrinsic x expected 
as a result of the emission of gravitational waves, 
x GR = —1.6 x 10~ 21 . Neglecting i int , the relation- 
ship between i and fi defined by equation || is pa- 
rameterized by the well determined physical param- 
eters x, x and fi. Also, because fi is fortuitously 
nearly anti-parallel to ft', 5i/5Q is close to zero, and 
incorporation of equation ^ in our timing model pro- 
vides a highly significant constraint on the inclina- 
tion angle. 

The orbital inclination parameterizes the shape 
of the Shapiro delay, that is, the delay due to the cur- 
vature of space-time about the companion. In highly 
inclined orbits, seen more edge-on from Earth, the 
companion passes closer to the line-of-sight between 
the pulsar and the observatory, and the effect is in- 
tensified. As the relative positions of the pulsar and 
companion change with binary phase, the Shapiro 
delay also varies and, in systems with small orbital 
eccentricity, is given by: 

A s = -2r ln[l - s cos(0 - <fo)]. (4) 

Here, s = sinz and r = Gm 2 /c 3 are the shape and 
range, respectively, cp is the orbital phase in radians, 
and (j) is the phase of superior conjunction, where 
the pulsar is on the opposite side of the companion 
from Earth (as shown in Fig. 1). For small inclina- 
tions, the orbit is seen more face-on from Earth, and 
As becomes nearly sinusoidal in form. 

In the PSR J0437-4715 system, the Shapiro ef- 
fect is six orders of magnitude smaller than the clas- 
sical Roemer delay, the time required for light to 




Figure 1 The three-dimensional orientation of the 
pulsar orbit is determined using a classical geometric 
model. With the centre of mass of the binary sys- 
tem at the origin, the basis vectors, lo, Jo and Ko, 
define east, north, and the line-of-sight from Earth, re- 
spectively. The orientation of the normal vector, n, 
is defined with respect to this basis by the longitude 
of the ascending node, fi, and the inclination angle, 
i. The plane of the sky, or lo-Jo plane, (shown stip- 
pled) intersects the orbital plane at the "line of nodes" 
(dashed line). Below the lo-Jo plane, the orbital path 
has been drawn with a dotted line. The unit vector, f2', 
lies in the lo-Jo plane and is perpendicular to the line 
of nodes. The pulsar is shown at superior conjunction, 
where radio pulses emitted toward Earth experience the 
greatest time delay due to the gravitating mass of the 
companion on the opposite side of the centre of mass. 

travel across the pulsar orbit. In nearly circular or- 
bits, the Roemer delay also varies sinusoidally with 
binary phase. Consequently, when modeling less 
inclined binary systems with small eccentricity, the 
Shapiro delay can be readily absorbed in the Roemer 
delay by variation of the classical orbital parameters, 
such as x. For this reason, a previous attempt at mea- 
suring the Shapiro effect in the PSR J1713 + 0747 
system 14 yielded only weak, one-sided limits on its 
shape and range. 

In contrast, we have significantly constrained the 
shape independently of general relativity, enabling 
calculation of the component of As that remains 
un-absorbed by the Roemer delay. The theoretical 
signature is plotted in Fig. 2 against post-fit residu- 
als obtained after fitting the arrival time data to a 
model that omits the Shapiro effect. To our knowl- 
edge, this verification of the predicted space-time 
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distortion near the companion is the first such con- 
firmation (outside our Solar System) in which the 
orbital inclination was determined independently of 
general relativity. 
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Figure 2 Arrival time residuals confirm the predicted 
space-time distortion induced by the pulsar compan- 
ion. The unabsorbed remnant of Shapiro Delay is 
much smaller than the theoretical total delay, which 
for PSR J0437— 4715 has a peak-to-peak amplitude of 
about 3.8/is. In the top panel, the solid line models 
the expected delay resulting from a companion with 
a mass of 0.236 M Q , at the geometrically-determined 
orbital inclination. Measured arrival time residuals, av- 
eraged in 40 binary phase bins and plotted with their 
la errors, clearly exhibit the predicted signature. In 
the bottom panel, the same residuals with the model 
removed have an r.m.s. residual of only 35 ns and a 
reduced x 2 °f 1-13. 

The range of the Shapiro delay provides an estimate 
of the companion mass, m 2 = 0.236 ± 0.017 M©, 
where M Q is the mass of the sun. Through the mass 
function 3 , f(M), we then obtain a measurement of 
the pulsar mass m p = 1.58 ± 0.18 M©. Slightly heav- 
ier than the proposed average neutron star mass 3 , 
TOp Vg = 1.35 ± 0.04 M Q , this value of m p suggests 
an evolutionary scenario that includes an extended 
period of mass and angular momentum transfer. 
Such accretion is believed to be necessary for a neu- 
tron star to attain a spin period of the order of a 
millisecond 15 . It is also expected that, during accre- 
tion, the pulsar spin and orbital angular momentum 
vectors are aligned. Under this assumption, the mea- 
sured inclination angle of i = 42? 75 ± 0?09 does not 
support the conjecture that pulsar radiation may be 
preferentially beamed in the equatorial plane 16 . 



The total system mass, M, can also be calculated 
from the observed Co, using the general relativistic 
prediction of the rate of orbital precession. Using 
M, f(M), and i, we obtain a second consistent esti- 
mate of the companion mass, m' 2 = 0.23 ± 0.14 M©, 
the precision of which is expected to increase with 
time as i 3 / 2 , surpassing that of the r-derived value 
in approximately 30 years. 

The complete list of physical parameters mod- 
elled in our analysis is included in Table 1. 
Most notably, the pulsar position, parallax dis- 
tance, d„ = 139 ± 3 pc, and proper motion, 
fj, = 140.892 ± 0.006 mas yr _1 , are known to accura- 
cies unsurpassed in astrometry. Although closer, d w 
lies within the 1.5 a error of an earlier measure- 
ment by Sandhu et al. 9 , 178±26pc. The d n and 
\x estimates can be used to calculate (3 and the in- 
trinsic spin period derivative, P mt = P ohs — (3P = 
(1.86 ± 0.08) x 10 - 20 , providing an improved char- 
acteristic age of the pulsar, r c = P/(2P mt ) = 4.9 Gyr. 
Another distance estimate may be calculated using 
the observed [i and Pb by solving Equation |3] for d, 
after noting the relative negligibility of any intrinsic 
contribution 17 . The precision of the derived value, 
c?b = 150 ± 9 pc, is anticipated to improve as t 5 ^ 2 , 
providing an independent distance estimate with rel- 
ative error of about 1% within the next three to four 
years. 

With a post-fit root mean square (r.m.s.) residual 
of merely 130 ns over 40 months, the accuracy of our 
analysis has enabled the detection of annual-orbital 
parallax. This has yielded a three-dimensional de- 
scription of a pulsar binary system and a new geo- 
metric verification of the general relativistic Shapiro 
delay. Only the Space Interferometry Mission (SIM) 
is expected to localize celestial objects with preci- 
sion similar to that obtained for PSR J0437-4715 
(including parallax) . By the time SIM is launched in 
2010, the precision of this pulsar's astrometric and 
orbital parameters will be vastly improved. Obser- 
vations of the companion of PSR J0437— 4715 using 
SIM will provide an independent validation and a 
tie between the SIM frame and the solar-system dy- 
namic reference frame. 

We also expect that continued observation and 
study of this pulsar will ultimately have an important 
impact in cosmology. Various statistical procedures 
have been applied to the unmodelled residuals of 
PSR B1855 + 09 (see ref. 19 and references therein) 
in an effort to place a rigorous upper limit on Q g , the 
fractional energy density per logarithmic frequency 



4 



W. van Straten et al. 



interval of the primordial gravitational wave back- 
ground. As the timing baseline for PSR J0437-4715 
increases, our experiment will probe more deeply 
into the low frequencies of the cosmic gravitational 
wave spectrum, where, owing to its steep power-law 
dependence 18 , the most stringent restriction on fl g 
can be made. 
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Table 1 PSR J0437-4715 physical parameters 

Right ascension, a (J2000) . . . 04 h 37 m 15!7865145(7 

Declination, S (J2000) -47°15'08"461584(8 

fi a (mas yr- 1 ) 121.438(6 

Us (mas yr" 1 ) -71.438(7 

Annual parallax, it (mas) 7.19(14 

Pulse period, P (ms) 5.757451831072007(8 

Reference epoch (MJD) 51194. ( 

Period derivative, P (10~ 20 ) . . 5.72 9 06(5 

Orbital period, P h (days) 5.741046(3 

x (s) 3.36669157(14 

Orbital eccentricity, e 0.000019186(5 

Epoch of periastron, T (MJD) 51194.6239(8 

Longitude of periastron, u> (°) . 1.20(5 

Longitude of ascension, £1 (°) . 238(4 

Orbital inclination, % (°) 42.75(9 

Companion mass, mi (Mq) . . . 0.236(17 

P b (10- 12 ) 3.64(20 

Co ( yr- x ) 0.016(10 



Best-fit physical parameters and their formal la er- 
rors were derived from arrival time data by minimizing 
an objective function, x 2 , as implemented in TEMPO 



( fittp: //pulsar. princeton.edu/tempo ). Our timing model is 
based on the relativistic binary model 19 and incorporates ad- 
ditional geometric constraints derived by Kopeikin 7,8 . In- 
dicative of the solution's validity, \ 2 was reduced by 30% 
with the addition of only one new parameter, Q. To de- 
termine the la confidence intervals of O and i, we mapped 
projections of the A% 2 = ^(Vl^i) — Xmin = 1 contour, 
where x 2 (f2, i) is the value of x 2 minimized by variation of 
the remaining model parameters, given constant SI and i. 
Parenthesized numbers represent uncertainty in the last dig- 
its quoted, and epochs are specified using the Modified Julian 
Day (MJD). 



